Order-disorder transition in nanoscopic semiconductor quantum rings 
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Using the path integral Monte Carlo technique we show that semiconductor quantum rings with 
up to six electrons exhibit a temperature, ring diameter, and particle number dependent transition 
between spin ordered and disordered Wigner crystals. Due to the small number of particles the 
transition extends over a broad temperature range and is clearly identifiable from the electron pair 
correlation functions. 

PACS numbers: 73.23.-b,73.20.Dx,31.15.Kb,71.45.Lr 



Nanoscopic semiconductor quantum rings (QRs), 
which recently have been experimentally realized by 
Lorke et. al Q], are next to quantum wires probably the 
best prototypes of quasi one dimensional quantum sys- 
tems. QRs can be viewed as rotating Wigner crystals 
with promising features for application in microelectron- 
ics. They can be modeled using a simple Hamiltonian of 
the form 
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where k = 12.9 and m* = 0.067me are the material con- 
stants of GaAs (|]||l2|. The parameter rp is the radius 
of the quantum ring and cjq defines the strength of the 
two-dimensional potential |^. Fig. ^(a) displays the 
shape of the ring potential. QRs can be tuned from quasi 
one-dimensional to two-dimensional systems by variation 
of the ring diameter and the potential strength. 




In this letter we present the results of path integral 
Monte Carlo (PIMC) simulations of single nanoscopic 
QRs with up to 8 electrons and different radii tq. We 
show that they undergo a temperature, radius and par- 
ticle number dependent spin order-disorder transition. 
Furthermore, the influence of quantum effects on the spa- 
tial electron distribution as well as the addition ener- 
gies Ai5 are given. Our results for the addition energies, 
i.e. the energy that is needed to place an additional elec- 
tron in a ring, is compared to the experimental results of 
Warburton et al. 0. 

In contrast to Hartree-Fock and spin density functional 
theory (DFT) PIMC samples without any approximation 
the full many body wavefunction instead of single or sums 
of Slater determinants. Especially for highly correlated 
electron systems this is a major advantage of PIMC. An- 
other benefit of PIMC is the possibility to study tem- 
perature dependent phenomena. For quantum dots the 
problems of different density functional approaches have 
intensively been discussed in the past |[ ||, The 
so called fermion sign problem is still a topic of actual 
research and restricts the application of PIMC to a lim- 
ited number of fcrmions and for QRs to a temperature 
of at least 10 K. 

The Feynman path integral for an A'^-electron sys- 
tem with position eigenket | r^, Si) [si — ±i for spin-up 
and spin-down electrons) in an external potential can be 
rewritten as llOtl 



FIG. 1: Potential energy surface of a semiconductor quantum 
ring (a) and schematic illustration of the definition of ip (b). 

While mesoscopic QRs have been investigated the- 
oretically and experimentally in depth Q, nanoscopic 
rings with strong quantum effects are of increasing in- 
terest. Koskinnen et al. Q reported CI calculations of 
rotational and vibrational excitations of nanoscopic QRs 
with up to N=7 electrons. They claim that QRs behave 
like rather rigid molecules or Wigner crystals with anti- 
ferromagnetic order in the ground state. Ahn et al. ||] 
considered stacked nanoscopic rings and found an N- de- 
pendent Stark effect. Experimentally it has been found 
that the emission energies of QRs change abruptly when- 
ever adding an electron 0. 
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and the boundary condition rj{M + 1) = rj{l). For 
M ^ cx) eq. (0) becomes exact. Standard Metropolis 
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Monte Carlo (MC) techniques can be utilized to evaluate 
the integral in (^). 

The basic quantities reflecting the spatial struc- 
ture of the electron configuration are the electron- 
electron (distance) pair correlation functions Tij{r) = 
{6{r—\rij I)), the angular pair correlation functions 
rij{(p) = {6{ip~\(pi — ipj\)), and the radial electron den- 
sity Pi{r) = 2^('^(''~ I from which all energies can 
be calculated using the Hypervirial theorem (for details 
of our method see |9| ) . The definitions of (p and are il- 
lustrated in Fig. ^b). To take the particle symmetry into 
account wc introduce Tij — F^-'- for Si ^ Sj, F^ = F^^^ 
for Si = Sj = i, and F^ = T^^ for Si — Sj — respec- 
tively. Obviously for S — we have T^^ — F^T 

In our simulations we controlled the systmatic error 
arising from the limited number of timeslices M and the 
statistical MC error carefully. By choosing M y.T ~ 600 
and using up to 10 biUion (10^°) MC steps per run we 
pushed the overall error of all energy expectation values 
below 0.3 %. Our Fortran-code is completely parallelized 
using MPI and Lapack and most calculations have been 
performed on a Cray T3E with 64 processors. 

We fixed the strength of the harmonic potential Tiloq — 
12 meV, resulting in effective atomic units for the Hartree 
energy E'^ = 10.995 meV and the Bohr radius Oq — 
10.1886 nm. A ring diameter rg — 14 nm — 1.37 Oq than 
corresponds to the experimental setup of Lorke et al.. To 
investigate the ring size dependence we performed addi- 
tional calculations for rg = 50.94 nm — 5aQ. 

Fig. ||(b) and (c) display the pair correlation func- 
tions of QRs with N — 6, S — and ring diameters 
To = 1.37 oji and 5.0 at T = 15 K. As expected for 
both diameters peaks at </? « 7r/3 and (f ~ 27r/3 occur 
in T{ip). For ro = 5ao the angular pair correlation func- 
tions for electrons with equal and unequal spin are almost 
identical, i.e. the Pauli principle does not play an impor- 
tant role in this case. In contrast, for tq ~ 1.37 Oq the 
pair correlation functions show a strong spin dependence. 
The electrons arrange on the ring with anti-ferromagnetic 
order. Such spin density waves have been predicted by 
Koskinen et al. |ll|. The role of quantum effects is 
reflected as well in Fig. (a). For the smaller ring size 



the radial electron density is much broader and nonvan- 
ishing at the ring center, implying that the system does 
not behave like a quasi ID-system. In a perfect hexago- 
nal Wigner crystal the equilibrium distances of the elec- 
trons would be r=l, 1.73, and 2.0ro. For ro — 5.0 Oq 
the distance pair correlation function is approximately 
a properly weighted superposition of Gaussians centered 
at these distances. For tq = 1.37 ag the pair correlation 
function depends on the total spin of a pair and is broad- 
ened up to ArQ. From Fig. ||(b) and (c) we infer that the 
most probable conflguration is one where the electrons 
are ordered on a zig-zag line around the circle, i.e. the 
electrons arrange alternately in the inner and the outer 
part of the ring. 

Fig. 1^ displays the angular spin density pair corre- 
lation functions for both ring diameters, temperatures 
T=15, 30, and 90 K and particle numbers iV=4, and 6. 
In all cases the large angle correlations disappear with 
increasing temperature and at T = 90 K only at small 
angles a spin correlation is still visible. For rg = 5.0 Oq 
the negative correlation at small angles increases with in- 
creasing temperature. However, this is simply due to the 
fact that the values of both correlation functions at small 
angles become larger with increasing temperature. The 
spin correlations for the small ring are about one order 
of magnitude larger than those in the ring with radius 
ro — 5.0aQ (Note the different scalings of the abscissae). 
A comparison between N — 4 and = 6 shows that the 
spin correlation is smaller for the larger system. A prob- 
able reason for this is that for A^ = 6 the contribution of 
the Coulomb repulsion is much larger (see Table ||) freez- 
ing the electron in the Wigner crystal and ~ thinking in 
the picture of one-particle wavefunctions - making the 
overlap between one-particle wavefunctions smaller. 

In summary, from Fig. |^ can be inferred that a spin 
order-disorder transition appears with increasing tem- 
perature, increasing electron number and increasing ring 
size. 

Next we consider the temperature dependence of the 
spin order-disorder transition in some more detail. Fig. ^ 
displays the mean values of the angular separation of the 
electrons with equal and unequal spin for A^ = 4 and 6, 
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FIG. 3: Angular spin density correlation function (r^^(<^) — 
rT^(yj)) for T = 10, 30, and 90 K and (a) = 4, 5 = 0, and 
ro = 1.37aS, (h) N = G, S = 0, and ro = 1.37aS, (c) iV = 4, 
5 = 0, and ro = 5a5, and (c) A^ = 6, 5 = 0, and ro = SaJ. 



5 = 0, and ro = 1.37 Oq as a function of temperature. 
As expected from the results presented above the over- 
all differences between the expectation values for equal 
and unequal spins approach zero with increasing temper- 
ature, i.e. the Pauli principle becomes less important. In 
addition, the values for = 6 are smaller than those 
for = 4. Obviously, this is because the available por- 
tion per particle of the ring volume is smaller for a larger 
number of electrons and the Coulomb repulsion is un- 
able to disperse the particles. At T = 10 K and = 4 
the contribution of the Coulomb term to the total poten- 
tial energy is 53 % for ro = 1.37 and only 36 % for 
ro = 5ao, while for = 6 the difference between the 
different diameters is with 60 % and 51 % substantially 
smaller (see Table |). 

The slope of {^p^-^) for A^=4 can be understood as fol- 
lows. Up to 40 K (v3Ti) 

increases due spin disordering. 
At higher temperatures {f^^) decreases due to increasing 
thermal fluctuations. 

Finally, we calculated the second energy differences 
IS.E = En+1 — 2En + En-1 also called addition ener- 
gies, which are an indicator of the stability of a quantum 
ring with a given number of electrons. For quantum dots 
it was claimed that the electron configurations are given 
by Hunds rule pH and consequently magic numbers oc- 
cur at A^ = 2, 6, 12, and 20. Here we calculated only 
the addition energy up to A/ = 7 for 25 K and N = 5 
for 10 K. As expected due to the Pauli principle and the 
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FIG. 4: Temperature dependence of the mean angles (0^^) 
and (0^^) for iV=4 and 6, S^O and ro = 1.37 aj. 



Wigner crystal structure a strong odd even effect occurs. 
The general behavior of this effect does not change for 
25 K. However, for higher temperatures significant differ- 
ences can be expected (see above). Warburton et al. 
argued that the general features of shell effects occur- 
ring in QRs are the same as in quantum dots. As can 
be inferred from Fig. |5| our calculations confirm this for 
N = 6. Furthermore, the addition energies are in the 
same range as those from photoluminescence measure- 
ments at 4 K 0. Due to the strong Coulomb repulsion 
in QRs, which grows with increasing electron number, it 
can be expected that shell effects become less important 
with increasing particle number. 




FIG. 5: Addition energies AE at temperatures T=10 and 
25 K for ro = 1.37 a^. 

The effects described above are reflected quantitatively 
in Table | presenting the total energy (-Etot), kinetic en- 
ergy (i?kin), total potential energy (i^pot), the energy due 
to the ring potential (-Ering), and the Coulomb energy 

In conclusion, we presented the results of full many 
body wavefunction calculations for QRs with up to eight 
electrons. We found that the properties of the rings de- 
pend in an intriguing manner on the ring diameter, the 
particle number and the temperature, which in turn is 
due to spin correlation. Coulomb ordering and the gen- 
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TABLE I; Total-, kinetic-, potential-, ring-, and Coulomb- 
energies for A'^=4 and 6, S = 0, tq — 1.37 and 5 for different 
temperatures. All energies are given in meV. 



N roK] 


r[K] 


-Etot 


-Ekin 


-Bpot 


-fi'ring 




4 1.37 


10 


57.2 


23.0 


34.2 


16.1 


18.1 




25 


59.3 


24.4 


35.0 


16.8 


18.2 




90 


85.1 


40.6 


44.4 


26.5 


17.9 


5 


10 


32.7 


15.0 


17.7 


11.4 


6.3 




25 


36.1 


17.7 


18.4 


11.6 


6.8 




90 


62.0 


35.4 


26.6 


18.5 


8.1 


6 1.37 


10 


117.8 


39.4 


78.4 


31.1 


47.3 




25 


122.5 


42.6 


79.9 


32.9 


46.9 




90 


160.9 


66.7 


94.2 


48.4 


45.8 


5 


10 


61.0 


25.3 


35.7 


17.3 


18.3 




25 


65.7 


28.8 


36.9 


17.8 


19.1 




90 


105.8 


55.3 


50.6 


29.1 


21.5 



eral strength of quantum effects. QRs exhibit a param- 
eter dependent spin order-disorder transition. By varia- 
tion of the ring diameter the system can be tuned from a 
quasi ID Wigner crystal to a 2D structure. The accesible 
parameter ranges can be used to tune the properties of 
quantum rings to desired values. Because of the ring di- 
ameter as an additional parameter this qualifies them as 
even better candidates than quantum dots for possible 
applications in microelectronics. The addition energies 
calculated using PIMC are in good agreement with the 
experimental results of Warburton et al. |7| and reflect 
the predicted shell effects. 
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